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Abstract 

Using the gauge/gravity duality we calculate the heavy quark-antiquark po- 
tential in a hot, anisotropic and strongly coupled Yang-Mills plasma in (3+1)- 
dimensions. As the anisotropic medium we take a deformed version of N = 4 super 
Yang-Mills theory at finite temperature following a recent work where the dual 
type IIB supergravity solution is also proposed. We turn on a small value of the 
anisotropy parameter, for which the gravity dual is known analytically (perturba- 
tively), and compute the velocity-dependent quark-antiquark interaction potential 
when the pair is moving through the plasma with a velocity v. By setting v = we 
recover the static quark-antiquark potential. We numerically study how the poten- 
tial is modified in the presence of anisotropy. We further show numerically how the 
quark-antiquark separation (both in the static and the velocity-dependent case) and 
hence, the screening length gets modified by anisotropy. We discuss various cases 
depending upon the direction of the dipole and the direction of its propagation and 
make a comparative study of these cases. We are also able to obtain an analytical 
expression for the screening length of the dipole moving in a hot, anisotropic plasma 
in a special case. 
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1 Introduction 

Ever since the seminal work of Maldacena [1], further elaborated in [2,3], the AdS/CFT 
correspondence (see [4] for a comprehensive review) and its subsequent generalizations 
have proved to be an indispensable tool for exploring the strongly coupled regime of large 
N gauge theories (where N is the rank of the gauge group). The correspondence, in its 
primitive incarnation, advocated the duality between type IIB string theory on AdS$ x S$ 
(where AdS stands for anti-de Sitter space) and J\f = 4, SU (N) super Yang-Mills (SYM) 
theory living on the 4-dimensional boundary of AdS^. Since then, it has been generalized 
to encompass a wider variety of gauge theories with known gravity duals and is now better 
called the gauge/gravity duality. One of the frontier areas where the duality has been 
particularly useful is the physics of quark-gluon plasma (QGP). 

Experiments at the Relativistic Heavy Ion Collider (RHIC) [5] have provided various 
fascinating insights into the properties of quantum chromodynamics (QCD) matter at 
extreme high temperatures, where it appears in the form of QGP. Most notably, the 
QGP does not behave as a weakly coupled gas consisting of quarks and gluons (as it 
should at very high temperature). Rather, there are strong indications that it resembles a 
strongly coupled fluid at the energy scale reached at RHIC [6] . This makes the theoretical 
description of QGP a challenging task. The strong coupling imposes severe limitations 
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on the applicability of standard perturbative QCD techniques. Lattice field theory is 
a powerful non-perturbative tool for exploring the static properties of strongly coupled 
gauge theories but has limitations when called upon to explain various real-time properties 
that are relevant to QGP. In such a scenario, the gauge/gravity duality has emerged as 
a promising tool for exploring such strongly coupled non-Abelian plasma ( [7] provides a 
review of the applications of gauge/gravity duality to QGP). Although the exact gravity 
dual to QCD has still eluded us, and the gauge theories and their gravity duals used for 
performing calculations are different from real world QCD, the results obtained so far have 
been quite encouraging. In fact, many of the results obtained exhibit a kind of universality 
among the different theories, the most notable among them being the celebrated rj/s 
ratio [8,9], where rj is the shear viscosity of the strongly coupled fluid and s is its entropy 
density. By now there is a large body of literature which attempts at calculating various 
quantities of interest, like the drag force, the jet quenching parameter, the static and 
the velocity-dependent quark- ant iquark potentials, the screening length, etc. in various 
QCD-like gauge theories in the deconfined phase using this duality. In many cases, the 
computed quantities have been in good qualitative agreement with their corresponding 
counterparts extracted from the experimental data. However, most of these works concern 
QGP which is locally isotropic. But, QGP, just after its creation in heavy ion collisions, is 
locally anisotropic and far away from equilibrium for a time t < r out . Further, it becomes 
locally isotropic only after time Ti SO > r out , so that the standard hydrodynamic description 
of the plasma makes sense only if we want to probe the plasma at time scale t > Ti so . 
To probe the QGP in the time window r out < t < r,, so , it is imperative that one takes 
into account the inherent anisotropy. It has been proposed [10-18] that an inherently 
anisotropic hydrodynamic description, which involves a derivative expansion around an 
anisotropic state, can be used to study the plasma in this regime. In this time domain the 
plasma has unequal pressures in the longitudinal and the transverse directions leading 
to an anisotropic expansion of the plasma. While in reality the degree of anisotropy 
will decrease with time, here we shall always treat that it is independent of time over 
a suitable time scale. Attempts to investigate the anisotropic plasma in the framework 
of the gauge/gravity duality include [19-40]. Recently, Mateos and Trancanelli [41,42] 
proposed a completely regular type IIB supergravity solution dual to anisotropic plasma. 
Further work involving this particular gravity dual can be found in [43-50] where the 
drag force, the jet quenching parameter, the stopping distance, the screening length, etc, 
were found out. In [51] it was shown that the 7]/s ratio violates the conjectured bound 
in anisotropic plasma. It is thus natural that one should further extend this program to 
explore the consequences of the presence of anisotropy and see how it affects the various 
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quantities of interest. Motivated by this, in this paper we address the issue of how the 
heavy quark-antiquark (Q-Q) potential is modified in the presence of anisotropy. We 
take the gravity dual as proposed in [41,42]. The static Q-Q potential and the screening 
length were already considered in [44] for such a system. Here we extend the analysis to 
the velocity-dependent case by considering a heavy Q-Q pair moving through the plasma 
with a velocity v. We consider only small values of the anisotropy parameter, in which 
case the dual gravity solution is known perturbatively. Following [52-55], we employ the 
gauge/gravity duality to first compute the expectation values of certain Wilson loops, 
which are nonperturbative objects in gauge theories. In the background given in [41], 
we introduce a fundamental string probe and extremize the Nambu-Goto string world- 
sheet action in a static gauge. This, in turn, yields the expectation value of the Wilson 
loop, where the loop is the boundary of the above minimal area. For the dipole velocity 
v < 1, the Wilson loop obtained is time-like and its expectation value can be related 
to the velocity-dependent interaction potential V of the dipole using the prescription 
proposed in [56]. We plot the potential V against the Q-Q separation L for various values 
of the velocity v and the anisotropy parameter a and study how the introduction of a 
small anisotropy influences the potential. Unlike in the static case (as described in [44]) 
where there were only two possible configurations of the dipole, here we shall see that the 
introduction of the velocity parameter gives rise to a plethora of possibilities, which we 
shall discuss one by one. We further probe the effect of anisotropy on the Q-Q separation 
and consequently the screening length L max . We are also able to obtain an analytic 
expression for the screening length in the anisotropic plasma in a special case. Although 
the static Q-Q potentials have already been provided in [44], we reproduce them here since 
they are recovered naturally in the v = limit of our analysis and nicely complements 
our results for the static Q-Q separation. In [49] the authors give an in-depth analysis 
of the screening length when the infinetely massive Q-Q pair (or the heavy meson or 
quarkonium, as they call it) moves in hot, anisotropic plasma. In this paper, one can 
read off the screening length from the plot of the Q-Q separation. Wherever the results 
overlap, they are in agreement with those obtained in [49]. 

The paper is organized as follows. In section 2, we briefly review the model and the dual 
geometry and discuss the general set-up. Section 3 is the core of the paper where we 
compute the Q-Q separation and the potential and provide the numerical results. We 
also calculate the screening length analytically in a special case. Section 3 is divided into 
five subsections corresponding to the different cases we consider. In section 4, we compare 
our results for the different cases considered and also with some other models available in 
the literature. Finally, in section 5 we summarize our work and conclude. 
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2 The dual geometry 



In this section we briefly review the gravity dual of the gauge theory we are interested 
in and discuss the general set-up of the problem. We will take the gauge theory as a 
deformed version of M — 4, SU(N) YM theory {N being the number of colors) at large 
t'Hooft coupling A = <?y M iV where the deformation is achieved by introducing a #-term 
in the action as 

S = Ssym + g^r J 0(x 3 )TtF A F (1) 

where 9(x 3 ) oc x 3 (we take (t . X - X - X ^ SiS the gauge theory coordinates). The presence 
of the non-zero #-term breaks the SO (3) rotational symmetry down to a 50(2) symmetry 
in the x 1 -x 2 plane and makes the theory anisotropic. The dual gravity background was 
given in [41,42] inspired by [57] and reads in the string frame, 

ds 2 = r 2 (^-FBdt 2 + (dx 1 ) 2 + (dx 2 ) 2 + H{dx 3 ) 2 + ^\ + e$+dnl (2) 

X = ax 3 , = 0(r) (3) 

where x is the axion, proportional to the anisotropic direction x 3 , the proportionality 
constant a being the anisotropy parameter and is the dilaton. r is the AdS radial 
coordinate with the boundary at r = oo and the horizon at r = r/j. Here dQ 2 is the 
metric on the five-sphere 5s and we have set the common radius of the AdS space and 
55 to unity. There is also a RR self-dual five-form which will not play any role in our 
discussion here. Anisotropy is introduced through the axion, the dual to the gauge theory 
#-term. The anisotropy parameter a turns out to be [41] a = Xti^/AwN where n^y is 
the density of D7-branes (which acts as the magnetic source of the axion) along the x 3 
direction. The D7-branes wrap around 5s and extend along the transverse directions, 
x l ,x 2 . Thus in the gravity dual the presence of anisotropy can be attributed to the 
existence of anisotropic extended objects. Note that the Z)7-branes do not extend along 
the radial direction. Hence they do not reach the boundary and do not contribute any 
new degrees of freedom to the theory. J 7 , £>, % are all functions of the radial coordinate 
r and are known analytically only in the limiting cases when the temperature is very 
high or low. Otherwise, they are known numerically in the intermediate range. J 7 is the 
'blackening factor' which vanishes at the horizon, i.e., jF{Th) = 0. The degree of anisotropy 
can be controlled by tuning the parameter a. In this paper, we shall be concerned with 
weakly anisotropic plasma (the small a or high temperature T limit, such that a/T <^ 1) 
in which case the functions J 7 , B, % can be expanded to leading order in a around the 
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black .D3-brane solution, 
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and we have defined the dimensionless quantity y = r jr^. The temperature is given by 

a 2 5 log 2 - 2 



T = — 



+ 0(a 4 



(6) 
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which can be inverted to yield the horizon position in terms of the temperature, which 
in the limit a/T 1, reads 

2 51og2-2 
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48tt 2 T 2 



(7) 



We introduce a fundamental string in this background and evaluate the Nambu-Goto 
action S. By extremizing this action we find the expectation value of the relevant Wilson 
loop. Assuming the string to move along x % with a velocity v and the string endpoints to 
lie along x\ separated by a distance L (which in the dual gauge theory translates to the 
quark- ant iquark separation), the Wilson loop so formed is a rectangle with a short side 
L along x- 7 and a long side T along any time-like direction in the t-x l plane. We further 
assume that L < T to ensure that the string world-sheet remains invariant under time 
translations. The expectation value of the Wilson loop so formed is, 



(W) 



AS -So) 



where Sq, the Nambu-Goto action corresponding to two disjoint strings (dual to a non- 
interacting quark and antiquark) is subtracted to offset the divergence in S. Now the 
Q-Q interaction potential V(L) is extracted from a knowledge of the expectation value of 
the Wilson loop via the working definition, 
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For the static Q-Q separation and potential one needs to consider only two possibilities: 
the dipole lying along the anisotropic direction x 3 or in the transverse plane. However, 
the introduction of the velocity v = tanh?7 (rj is the rapidity parameter) opens up the 
following possibilities: 

1. Motion in transverse plane, dipole lies perpendicular to direction of motion in the 
plane. 

2. Motion in transverse plane, dipole along x 3 . 

3. Motion along x 3 , dipole in transverse plane. 

4. Motion in transverse plane, dipole parallel to direction of motion in the plane. 

5. Motion and dipole, both along the anisotropic direction. 3 

3 Q-Q separation and Q-Q potential 

In this section we discuss the different cases, alluded to above, one in each subsection, and 
for each case we numerically study the Q-Q separation with varying values of the rapidity 
parameter rj and the anisotropy parameter a and see how the separation and hence, the 
screening length gets affected when we turn on a small value of a. We also compute the 
Q-Q potential (both velocity-dependent and static) and observe the modifications brought 
about by anisotropy. Further, we provide an analytic expression for the screening length 
in a special case. 

3.1 Motion in transverse plane, dipole lies perpendicular to di- 
rection of motion in the plane 

In this case the motion is wholly contained in the transverse plane and the dipole presents 
itself perpendicular to the direction of its motion. We first set our axes such that the dipole 
moves along x 1 while itself being aligned along x 2 . Then we go to the rest-frame (t',x v ) 
of the Q-Q pair via the following coordinate transformation, 

dt = cosh t)dt' — sinh rjdx 1 ', 

i ( 10 ) 
dx 1 = — smhrjdt' + coshrjdx ' . 

3 Of course, there exist other possibilities where the dipole can have any arbitrary orientation with 
respect to its direction of motion, which, itself, can be in any arbitrary direction. However, we do not 
consider these cases here. 
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Now the Q-Q pair and hence the Wilson loop can be regarded as static in a plasma that 
is moving with a velocity v in the negative x v direction. This implies that the rectangular 
Wilson loop spans the t' (since x v is fixed in this rest-frame) and x 2 directions with sides 
T and L respectively. In terms of the boosted coordinates the metric (2) can be rewritten 

as 



ds 2 = -A{r)dt 2 - 2B(r)dtdx 1 + C(r)(dx 1 ) 2 + r 2 ( [dx 2 ) 2 + U(dx 3 ) 2 + 
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(Note that since we shall be using the primed coordinates from now on, we have got rid 
of the primes for simplicity. Also, we have suppressed the ^-dependence of the quantities 
T% £>2- Further, we have expressed A, B and C as functions of the scaled radial coordinate 
y.) In this background we evaluate the Nambu-Goto string world-sheet action, 
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where g a p is the induced metric on the string world-sheet, 
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Here is the ten-dimensional metric as given in (11) and £ Q,/3 are the world-sheet 
coordinates, £° = r and = a. We choose the static gauge for evaluating (13) as 
r = t, a = x 2 where —L/2 < x 2 < +L/2 and r = r(a),x 1 (a) = x 3 (a) = constant. We 
wish to determine the string embedding r(a) supplemented by the boundary condition 
r (x 2 = ±L/2) — > oo. Equipped with the above parametrization, the Nambu-Goto action 
(13) can be written as, 
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At this stage, for convenience, let us define the new dimensionless quantities, a = ajr^ 
and I = L/rh, the scaled Q-Q separation, whence we can rewrite the above action as 

1/2 

S = Ilh. f ddC (16) 
2vra' J y ' 

-1/2 

where 



C=y/A(G22 + G rr y*) (17) 

is the Lagrangian density and d&y = y'. Note that C does not have any explicit in- 
dependence which at once allows us to extract the conserved quantity, 

,dC AG 22 T ^ , , 

£ - y — = . = K 18 
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which, in turn, yields, 
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Upon integration we obtain 

1/2 00 

y t 

The limits in the second integration require a little explanation. Recall that y is the scaled 
radial coordinate and the string hangs down starting from y = 00 (where the boundary 
gauge theory lives) up to yt (which we shall find shortly), where it turns back and rises 
again up to y = 00. Plugging in the explicit expressions for the metric components, we 
arrive at, 

~ 00 

2,k r 1 

l = — dy (21) 
Th 1 V (2/ 4 - 1 + tj£(y)) (V 4 - yt + |A(y) cosh 2 rj) 

where we have defined, 

S(y) = 8(y 2 - 1) - 10 log 2 + (3y 4 + 7) log M + lj , (22) 

A(y) = 2(1 - y 2 ) - 10 log 2 + 2(y 4 + 4) log + ^ J , (23) 
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and a = a/r^(~ a/7rT),K 
actual Q-Q separation as 



Kjr\,y\ = cosh 2 77 + K 2 . Using (7) one can now find the 




00 
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(24) 



As mentioned earlier, to perform the integration, one needs to specify yt. The turning 
point is found out by demanding that the terms in the denominator vanish separately 
(which is equivalent to demanding that y' vanishes at these points) and accepting the 
larger one among them. As one can easily verify, the first term in the denominator 
vanishes at y = 1, since = 0, thereby, furnishing a turning point at y t \ = 1 up 

to 0(a 2 ). To find the turning point y t 2 arising from the second term, we assume the 
anisotropy parameter a 4 to be small and we need to find a solution to 



Note that we have evaluated A at y = y c since the term is already at 0(a 2 ) and conse- 
quently, the error incurred is ~ (9(a 4 ). This has a solution 



It can be shown that yt2 > 1 always, so that we take it to be the actual turning point yt. 
As expected, by setting a = we recover the turning point y c in the isotropic case. Now 
(24) gives the Q-Q separation as a function of the constant K. However, it is not possible 
to perform the integration analytically. Hence, we shall resort to numerical integration 
and show how the Q-Q separation is affected by the presence of anisotropy for various 
values of the rapidity parameter 77 and the anisotropy parameter a. Further, by setting 
77 = we obtain the static Q-Q separation. We postpone the discussion of our numerical 
results till we give the Q-Q potential. 

Changing the integration variable from a to y we can rewrite the action (16) as 
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= a 2 , we shall, henceforth, call a the 
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Putting the explicit expression for the metric components one finally has 
Tr h f , y 4 - cosh 77 + |jA(y) cosh rj 



S = — ? / dy 

™ i v / (z/ 4 -l + Ss( 2/ ))(^- 2/ 4 + |A( z/)cos h 2 r / ) 

00 

= 7J± f dyS ani . (28) 

wt 

A mere inspection reveals that the integral, as it stands, diverges. This is because the 
action 5" contains the self-energy contributions from the free heavy Q-Q pair which, them- 
selves, are divergent. To obtain the Q-Q interaction potential V(L) we need to cure this 
divergence, which is done by subtracting from S the action So of a free Q-Q pair whence 
from (8,9) 

V(L) = (29) 

To compute So we consider an open string hanging down the radial direction (in gauge 
theory it translates to a single quark/ant iquark propagating in the same background as 
before) and employ the static gauge condition, r = t, a = r, x l = x l (a) and x 2 ,x 3 are 
independent of r, a. We evaluate the Nambu-Goto action (and multiply by two to take 
into account the contribution from the quark and the antiquark) which takes the form, 

00 

S = — t I dr^/AGrr + ix^iAC + B 2 ) 
iia' J 
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As before, So too does not have any explicit x 1 -dependence implying that there exists a 
conserved quantity, 
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(32) 



(we have used the scaling Kq = Ko/r 2 ). Note that in the expression for So we have 
not specified the lower limit of the integration ro (or y$ after scaling), which we shall 
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now determine. For a string (corresponding to a free quark/ant iquark) hanging down 
we expect it to extend all the way to the horizon at y = 1. This is the case when the 
string moves through the isotropic background. In particular, this implies that the string 
can not encounter a turning point before y — 1. In our case, the possible turning points 
can be found out from (32) by demanding that x 1 ' = oo at those points. Now the first 
two terms in the denominator of (32) give the turning point yo = 1 up to 0(d 2 ) since 
S(l) = A(l) = 0. However, the third term (which contains the unspecified constant K ) 
gives a turning point y$ ~ 1 + Kq + 0(d 2 ) which is greater than zero even for the isotropic 
case. Taking cue from the isotropic case we eliminate this possibility by constraining the 
value of Kq such that the zero of this term coincides with the zero of the numerator. This 
at once provides us an expression for Kq as 



K 2 = sinh 2 7] ^1 - (y = ^coshr/)^ . 



(33) 



Thus we conclude that even in the presence of anisotropy the string extends right up to 
the horizon (like in the isotropic case) without picking up any correction at least to 0(d 2 ). 
We can now recast the action as 
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Inserting (28) and (34) in (29) and then using (7) we can now write 

(oo 
I dyS° 




(35) 



where we have used the standard AdS/CFT dictionary R A = \a' 2 (with R set to unity 
here) to express our final result in terms of quantities pertaining to the gauge theory. 
Evaluating (35) involves performing integrals which can not be handled analytically. We, 
therefore, fall back upon numerical means to perform these integrals and numerically 
show our results. We compute the Q-Q separation L for various values of the rapidity 
parameter rj and the anisotropy parameter a as a function of the constant K, numerically 
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Figure 1: (a) shows the plot of quark-antiquark separation L (normalized) as a function of 
integration const. K with r/ = 1 for different values of the anisotropy parameter a when the 
dipole lies perpendicular to its velocity but both lie in the transverse plane, (b) shows the plot 
of properly normalized quark-antiquark potential V as a function of L with rj = 1 for the same 
set of anisotropy parameter values. 



invert (24) to express K in terms of L and plug it in (28) to finally obtain the Q-Q 
potential as a function of the Q-Q separation. Here we shall provide our numerical results 
for both the Q-Q separation and the Q-Q potential. In particular, by setting rj = we 
recover the static Q-Q potential, which was already found out in [44]. In Figs.l and 2 we 
have presented the Q-Q separation L(K) as a function of K and the Q-Q potential V(L) 
as a function of L for rj — 1 (v — 0.762) and r\ = 4 (v = 0.999) respectively. Each figure 
has two parts - part (a) showing the L(K)-K plot and part (b) showing the V(L)-L plot. 
While the qualitative pattern of the plots in both the figures are the same, the details 
differ. So we shall take Fig.l as the prototype case and discuss the results. First of all, 
we find from part (a) that as K increases the separation L(K) increases till it reaches a 
maximum L max after which it again falls off. L max is interpreted as the screening length 5 
of the dipole, i.e., beyond this critical value of L the screening effect of the plasma is 
sufficient to break the dipole. We observe that the effect of anisotropy is to suppress the 
screening length thereby encouraging the melting of the dipole. In particular, the degree 

5 Note that our definition of the screening length differs slightly from that used in [49]. 
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Figure 2: (a) shows the plot of L (normalized) as a function of K with r\ = 4 for different values 
of a when the dipole lies perpendicular to its velocity but both lie in the transverse plane, (b) 
shows the plot of properly normalized V as a function of L with r\ = 4 for the same set of values 
of a. 

of suppression of L max is more for stronger anisotropy. The deviation from the isotropic 
curve is more pronounced for lower K (before L max is attained) than for higher K (after 
Lmax)- For L < L max there can be two dipoles at a fixed L for two different values of 
K. To understand at which one of the K values the dipole will actually exist we need 
to analyze the V(L)-L plot. The Q-Q potential has two branches corresponding to two 
different values of K. The upper branch corresponds to smaller value of K whereas the 
lower branch corresponds to higher value of K. Of course, the lower branch has lower 
energy and consequently, it is the preferred state of the dipole. So, even if a dipole is 
in the upper branch it will not be in a stable configuration and the dipole will make a 
transition to the lower branch. As we turn on a small anisotropy both the branches of 
the potential shift slightly upwards. Since the upper branch is physically insignificant, 
corresponding to an unstable state, we shall confine our discussion to the lower branch 
only. The marginal upward shift in the potential indicates that the dipole is now loosely 
bound, the shift being more prominent for higher values of the anisotropy parameter 
a. This fits in with our conclusion from the L(K)-K plot that anisotropy enhances the 
screening effect of the medium. Further notice that in both cases the potential is always 
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Figure 3: (a) shows the plot of L (normalized) as a function of K with r\ = for different values 
of a when the dipole lies in the transverse plane, (b) shows the plot of properly normalized V 
as a function of L with 7] = for the same set of values of a. 



negative. In Fig. 3 the plots for the static case are shown. Since the basic nature is the 
same, we shall not elaborate upon our results and briefly mention the salient features 
of the plots emphasizing the differences from the velocity-dependent cases. First of all, 
notice that unlike the moving dipole case, now the deviation from the isotropic curve in 
the L(K)-K plot is appreciable on either side of L max . Also note that now L max is much 
higher for the rj = case and steadily decreases as we increase a. In the V(L)-L plot the 
lower branch suffers a small elevation whereas, the insignificant upper branch is largely 
insensitive to changes in a. The new feature that now emerges is that the static potential 
crosses zero and becomes positive at a particular value L = L p . As is well-known, in the 
presence of the medium the potential has two parts - the Coulomb part (varying inversely 
with L) and the confining part (which goes as L). L p denotes the separation beyond 
which the confining part starts to dominate. This feature is nicely captured in the plot 
here in a qualitative manner. In fact, there will be a critical velocity v p = tanh?7p (whose 
value will, in general, also depend upon a) beyond which the potential will not contain 
any positive piece. 

We shall now obtain an analytical expression for the screening length, albeit in a special 
case. The Q-Q separation L has already been given in (24). We have mentioned earlier 
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that, in general, the integration appearing in (24) can not be done analytically. Of course, 
this is not to be thought of as the artifact of our anisotropic background. Rather, it is a 
handicap present in the isotropic case, too. Here, to facilitate analytical manipulation, we 
shall confine ourselves to the ultra-relativistic regime where 77 is large, in which case the 
turning point y t also becomes very large but assume the product a 2 cosh 2 77 is sufficiently 
small. In this special case the first term in the denominator in (24) lends itself to a 
binomial expansion. Here, for the sake of simplicity, we shall consider only the leading 
order term in the aforesaid expansion in which case one can write 



hi 



dy = +... (36) 



Also, in the limit 77 becoming very large, A(y) reduces to A(y) = 1 — 10 log 2, which is, in 
fact, independent of y. In this simplified scenario, the integral can be handled analytically 
and we have, 

r 2k A . 5 2 , K ,_o v^r(3/4) 



-(5 log 2 -2) Ph; — (37) 
ttT V 48 1 ') yl T(l/4) 1 ; 

It is now a straightforward exercise to compute the value of K and hence, yt which 
maximize L as 

~ 2 n l2 a 2 cosh 2 r] (10 log 2 - 1) 
K = 2 cosh 2 r/ + 



Vt - ^ 1+ Wg2-l) cos hS). (3. 

Incorporating these values we arrive at the final expression for the screening length L mc 

as, 

1 1X3/4)2^ 1 / a 2 , l9 A 

v^FT T(l/4) 3 3 / 4 v ' \ lQ\(l-v 2 ) 
The proportional change brought by anisotropy is, 



(39) 



AL max ~a 2 , 2 

- — = 2.965735 cosh 2 n - 0.488578 . 40 

L I 16 

max 

I a=o 

Having deduced the analytical expression for L max , a few comments are in order here. 
First, as expected, by setting a = here one recovers the usual screening length in an 
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Figure 4: (a) shows the plot of L (normalized) as a function of K with 77 = 1 for different values 
of a when the velocity is in the transverse plane and dipole lies along anisotropic direction, (b) 
shows the plot of properly normalized V as a function of L with 77 = 1 for the same set of values 
of a. 

isotropic plasma [7,58]. Second, it is obvious that the correction factor is always negative 
so that L max decreases in the presence of anisotropy. Third, when a increases, the fall 
in L max is greater. Again, keeping a fixed, if r\ increases, L max falls. These conclusions 
drawn from the analytic expression (39) are in agreement with all our numerical results in 
the L(K)-K plots discussed earlier. Observe that the correction in the screening length 
arising due to the presence of anisotropy depends on the rapidity parameter as well. One 
also finds that L max depends inversely upon the temperature and scales with velocity as 
(1 — w 2 ) 1 / 4 . The velocity-scaling obtained here is in agreement with that found in [49], 
where, of course, arbitrary orientation of the dipole with respect to its velocity was allowed 
and the analysis was not restricted only to weak anisotropy. One infers from (38) that 
the value of K which maximizes L increases when we turn on the anisotropy parameter. 
This is also nicely exposed in the L(K)-K plot in Figs.l and 2 where the peaks gradually 
shift towards right as the anisotropy gets larger. With this we close our discussion of this 
configuration and move over to the next case. 
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Figure 5: (a) shows the plot of L (normalized) as a function of K with r\ = 4 for different values 
of a when the velocity is in the transverse plane and dipole lies along anisotropic direction, (b) 
shows the plot of properly normalized V as a function of L with r\ = 4 for the same set of values 
of a. 



3.2 Motion in transverse plane, dipole along x 3 

In this case the dipole lies along the anisotropic direction x 3 and moves in the transverse 
plane with a velocity v. Without any loss of generality, we can take the direction of 
motion to be along x 1 . The calculation in this case proceeds in pretty much the same 
way. So we shall be brief in this section, pointing out only the differences that crops up 
in the calculations as we go along. Firstly, note that the choice of the static gauge is 
slightly altered. Now we take r = t,a = x 3 ,r = r(a) with x 1,2 being independent of r or 
a. Enforcing this choice of gauge in the Nambu-Goto action (13) results in the following 
form of the action 

+1/2 

S = ~^t< S dd ^^ A ^ + G "-y' 2 ))- ( 41 ) 

-1/2 

As before, the absence of any explicit a-dependence leads to the following conserved 
quantity, 

K = AGs3 (42) 

^A(G 33 + G rr y' 2 ) 
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(a) (b) 

Figure 6: (a) shows the plot of L (normalized) as a function of K with 77 = for different values 
of a when the dipole lies along anisotropic direction, (b) shows the plot of properly normalized 
V as a function of L with rj = for the same set of values of a. 



and the scaled Q-Q separation assumes the form, 



, = !ur /*feas5hfs- (43) 

yt 

Plugging in the explicit expressions of the metric components, we finally obtain the Q-Q 
separation as 

_ 00 



nT V 48 J t \j{y' - 1 + £e) (y 4 - (i - I A) cosh 2 v - km-' 

(44) 

where we have suppressed the explicit ^-dependence of £(?/), A(y) and T-L(y). The turning 
point yt is found out by demanding the vanishing of the second term in the denominator 
at y t , i.e., from 

a 2 

yt - cosh 2 n + T^Myt) cosh 2 7i - K 2 U{ytT x = 0. (45) 
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The action is given by 



s -S/WlCT (46) 



Hi 



which can be rewritten as 



Tr h f y 4 - cosh 2 77 + §j A cosh 2 77 



Tr h 



(y*-l + gS) (y 4 - (1 - |A) cosh 2 77 - ) 
d7/S an \ (47) 



.y/ 



To evaluate the Q-Q potential one also needs to subtract the self-energy term S - It is 
easy to convince oneself that in this case the expression for £0 as given in (34) remains 
unaltered and the Q-Q potential will be given by (35) with S am now taken to be as in 
(47). We have given the L(K)-K and the V(L)-L plots in Figs. 4 and 5 for rj = 1 and 
77 = 4 respectively. Fig. 6 shows the static Q-Q separation and the static Q-Q potential. 
We observe that in all the cases the general pattern of the plots (like the rightwards shift 
of the peak in the L(K) curves, attenuation of L max and rise in the V(L) plots with 
increasing a) mimic those obtained earlier in Sec. 3.1 and hence does not merit a separate 
discussion. 

3.3 Motion along x 3 , dipole in transverse plane 

Third in our list is the case where the dipole is aligned in the transverse plane and it has 
a velocity along the anisotropic direction. For the sake of simplicity we have taken the 
dipole to lie along x 1 . While we shall proceed along the same line as in the previous cases, 
this time the calculations will be a little different since we now need to give a boost along 
the anisotropic direction, a; 3 . First of all, we go to the rest-frame (t',x 3 ') of the Q-Q pair 
by inflicting the boost 

dt = cosh rjdt' — sinh ridx 3 ' 

dx 3 = — sinh rjdt' + cosine/a; 3 '. 

The Wilson loop so formed spans the t' and x 1 directions. In terms of the boosted 
coordinates the metric (2) can be rewritten as 

ds 2 = -A(r)dt 2 - 2B(r)dtdx 3 + C(r)(dx 3 ) 2 + r 2 ((dx 1 ) 2 + (dx 2 ) 2 + + e^rffi 2 

= G^dx^dx" (49) 
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Figure 7: (a) shows the plot of L (normalized) as a function of K with rj = 1 for different values 
of a when the velocity is along anisotropic direction and dipole lies in the transverse plane, (b) 
shows the plot of properly normalized V as a function of L with 77 = 1 for the same set of values 
of a. 



where 



My) 

B(y) 
C(y) 



a 2 



y A — cosh 2 rj + —A cosh 2 rj + y A sinh 2 77(1 — H) 



sinh rj cosh rj 



a 2 



i-^A + y 4 (H-i; 



a 2 



y 4 + sinh 2 rj — —A sinh 2 rj + y A cosh 2 T]^H — 1) 



(50) 
(51) 
(52) 



To evaluate the Nambu-Goto string world-sheet action we employ the following choice of 
gauge: r = t,a = x 1 ,r = r(a) with x 2,3 having no r- or a- dependence. The action (13) 
can now be written as 

+1/2 



S 



Tr h 
2na' 



da \ A Gu + GrrU 
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(53) 



-1/2 



Again the absence of any explicit a-dependence furnishes the conserved quantity, 



(54) 



A Gu + G rr y 
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Proceeding in the same way as before we get the scaled Q-Q separation, 



,=at '<^ (55) 



m 



from which one can read off the actual Q-Q separation 



= ■ (56) 
V - Vt + fi A cosh2 ^ + 2/ 4 sinh2 ^{l-H) 

The turning point |/t is found from the solution of 

Vi ~ vt + ^Myt) cosh 2 rj + y 4 sinh 2 77(1 - H) = 0. (57) 

The factor (1 — "H) goes as ^ log ^1 + ^\ up to 0{a 2 ) and for large ?/ its contribution 

to the second factor in the denominator is ^-y sinh 77. This is greater than the other 
anisotropic term by 0(y 2 ) for large y. Since, the integration over y extends up to y = 00, 
the contribution from this term can be quite large. Hence, unlike in the previous cases, 
this time we do not expect the turning point yt to appear in the form of a correction 
to the isotropic value y c since the presence of this 0(y 2 ) term renders the applicability 
of perturbative methods to solve the above equation futile. Thus, one has to depend 
solely upon numerical techniques to solve (57) in order to extract y t . In fact, numerical 
evaluation shows y t to be markedly different from y c , particularly for low values of K. 
Once we have obtained y t , we use it in (56) to numerically study the Q-Q separation. 
The string world-sheet action is 



dyAj °" g " (H 

7T«' J V AG U - K 



Vi 



which, written explicitly, assumes the following form, 

Tr h f 1 y 4 — cosh 2 77 + cosh 2 77 + y 4 sinh 2 7/(1 — T-L) 



S = — ? dy 



na 

yt 
00 



(y 4 - 1 + fjE) (t/ 4 - yt + |A cosh 2 7/ + t/ 4 sinh 2 r/(l - H)) 
dj/5 ani . (59) 



m 
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Figure 8: (a) shows the plot of L (normalized) as a function of K with 77 = 3 for different values 
of a when the velocity is along anisotropic direction and dipole lies in the transverse plane, (b) 
shows the plot of properly normalized V as a function of L with 77 = 3 for the same set of values 
of a. 



As in the preceding cases, this action is divergent which is cured by taking away the 
self-energy contribution Sq of the Q-Q pair. To compute So we consider an open string 
hanging down the radial direction in the following gauge, r = t,a = r,x 3 = x 3 (a) and 
x x ,x 2 are independent of r, a. Repeating the same exercise as in Sec. 3.1 one finds Sq to 
be of the form 



T f J AC + B 

So = — : / dr\ AG 



ira' J v \{AC + B 2 -Kl 



Trh °f v/y 4 — cosh 2 77 + cosh 2 77 + ?/ 4 sinh 2 77(1 — "H) 

— 7 / d V / — x 

™ { Jy^-l-kl + ^KU+iH-lW-l) 



yt-\ + *m + {u-iw 



Tr h 



dySt (60) 
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where Kq is the conserved quantity owing its origin to the absence of any explicit in- 
dependence in the action. The second terms each in the numerator and the denominator 
separately vanish at y — 1 providing a potential turning point yt — 1. The first term 
in the denominator can contribute another turning point yt > 1 but that possibility is 
ruled out by judiciously choosing the constant K such that the zero of the first term in 
the numerator coincides with that of the first term in the denominator. We are now in 
a position to finally compute the Q-Q potential (35) with S an% provided in (59) and the 
corresponding self-energy term S^ 11 in (60). Using the above information we have plotted 
the Q-Q separation and the Q-Q potential in Figs. 7 and 8 6 . While the gross features 
of the plots remain almost unaltered, observe that all the signatures of the presence of 
anisotropy are far more pronounced (particularly in the high rapidity regime) than in 
either of the preceding cases. This has its roots in the presence of the 0(y 2 ) term in 
the anisotropic contribution to the Q-Q separation and the Q-Q potential as mentioned 
earlier. The heavy Q-Q potential for this configuration has also been found in [49], using 
different values of the parameters and we find that our results tally with those presented 
in [49]. 



3.4 Motion in transverse plane, dipole parallel to direction of 
motion in the plane 

We now come to the case where the dipole is aligned parallel to its direction of motion. 
This common direction can be in the transverse plane or along the anisotropic direction. 
We consider the former case in this subsection. For simplicity we shall take this common 
direction to be along x 1 . Boosting to the rest-frame, choosing the static gauge, r = t,a = 
x 1 , r = r(a) and x 2 = x 3 = constant leads us to the action, 

1/2 

S = S^ J ^VAiC + Grry^ + B 2 (61) 

-1/2 

which, in turn, supplies the constant of motion, 

AC + B 2 

K = . 62 
^A(C + G rr y' 2 ) + B 2 



6 Note that we have not given the static Q-Q separation and the static Q-Q potential in this case since 
these will be the same as in Sec.3.1. 
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Figure 9: (a) shows the plot of L (normalized) as a function of K with rj = 1 for different values 
of a when the dipole is parallel to its velocity and both lie in the transverse plane, (b) shows 
the plot of properly normalized V as a function of L with r/ = 1 for the same set of values of a. 



Proceeding along the lines of the earlier cases, we compute, 



y'-i 1 / , 

K Jy 4 - cosh 2 ?] + f^Acosh 2 ?] 



vrT V 48 



from which we find the Q-Q separation to be 

2K ( 5 2 (51og2-2)\ 7 , V - cosh 2 + |A cosh 2 r/ 

1 + - a / d v i x 

J l yV - 1 + W) (v 4 - 1 + |a) 

i 



y 4 -i-^ 2 + |A 



24 

The turning point y t is obtained from (63) which satisfies 



(64) 



yf-l-K 2 + ^A (y t ) = 0. (65) 

At the same time, note that y' now encounters a singularity at y s , given by, 

a 2 

y A s — cosh 2 r/ + — A(y s ) cosh 2 77 = 0. (66) 
25 




Figure 10: (a) shows the plot of L (normalized) as a function of K with rj = 4 for different 
values of a when the dipole is parallel to its velocity and both lie in the transverse plane, (b) 
shows the plot of properly normalized V as a function of L with 77 = 4 for the same set of values 
of a. 



Further, it is evident that for y < y s , the numerator in (64) becomes imaginary. So any 
potential turning point has to satisfy 

a 2 

y\ — cosh 2 7] + —A(y t ) cosh 2 rj > (67) 
which imposes a lower bound on K that turns out to be 7 , 

K 2 > K 2 mm = sinh 2 n - (y = Vcoshr/) ) . (68) 

Incidentally, note that this lower bound turns out to be the same as the constant Kq that 
appeared in Sec.3.1. Upon simplification the action boils down to, 

The existence of this lower bound is found in the isotropic case too as given in [56] . 
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oo 

S=TTh f dy^AG r 



AC + B 2 
AC + B 2 - K 2 



Tr h 7 vV " ^sh 2 rj + S A cosh 2 rj) (y 4 - 1 + f A) 
dy— 



yt 



Tr h 



(2/ 4 -i + |£) (y 4 -i-^ 2 + gA 

dyS am (69) 



with K respecting the inequality in (68). The self-energy contribution Sq is also given 
by (69) but now with K saturating the bound in (68) so that Sq becomes identical with 
that given in (34). We can now compute the Q-Q potential using (34) and (69) in (35) 
with (7). The Q-Q separation and the potential have been plotted in Figs. 9 and 10 for 
rj = 1 and rj = 4 respectively. The L(K)-K plots show that curves for higher value of 
K (after L max is attained) exhibits the same pattern as in the earlier cases but for lower 
values of K there is an inaccessible region for K < K m i n for which there is no solution 
to the dipole separation. This is reflected in the V(L)-L plot where the upper branch of 
the potential terminates abruptly at L — L min whereas the lower branch shows the usual 
behaviour . A closer scrutiny of the figures suggest that K min increases with increasing 
a and concomitantly, L m i n decreases. However, this is manifested only in the unstable, 
high energy branch, which, in any case, is devoid of much physical significance, being 
energetically unfavorable. 

3.5 Motion and dipole, both along the anisotropic direction 

Finally, we take up the case where the dipole is oriented along the anisotropic direction 
x 3 and it moves in the same direction. This time we shall make use of the metric (49) as 
obtained in Sec. 3. 3 and use the gauge choice of Sec. 3. 2. All the calculations proceed in 
identically the same fashion as in Sec. 3.4 and we end up with the Q-Q separation 



2K ( 5 2 (51og2-2)\ f \/2/ 4 -cosh ry + f^Acosh r/ + y 4 sinh r]{l-U) 

1 + jz / d v- — = — x 



7rTV 48 JJ yt yV-l-i^ + SA^ + ^-lX^-l) 

1 = (70) 
(y 4 - 1 + g£) (y 4 - 1 + |AH + (H- l)(y 4 - 1)) 

The first term in the denominator provides the turning point y t which, in turn, is con- 
strained by the condition that the numerator must be real. This results in a lower cut-off 
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(a) (b) 

Figure 11: (a) shows the plot of L (normalized) as a function of K with 77 = 1 for different 
values of a when the dipole is parallel to its velocity and both lie along the anisotropic direction, 
(b) shows the plot of properly normalized V as a function of L with rj = 1 for the same set of 
values of a. 



on the value of K. We can read off this lower bound by demanding that the zeros of the 
numerator and the first factor in the denominator occur at the same value of y. As was 
the case in Sec. 3. 3 due to the presence of the y 4 sinh 2 r](l — 7i) term here we do not expect 
the effect of anisotropy to be small enough so as to employ perturbative methods. Hence, 
we have evaluated the lower limit K min and y t completely numerically. Finally, the action 
becomes, 



Tr h f yV 4 ~ cosh rj + ^Acosh 77 + |/ 4 smh 77(1 — %) 
S = / dy- 



Ka ' I Jy*-l-ki + ^kH + {H-lW-l) 



y 4 -l + SAH + (H-l)(y 4 -l) 



/ 4 - 1 + §S 



21" 



not' 

yt 



dyS ani . (71) 



In similar fashion, one finds the self-energy contribution So to be the same as in (60) 
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Figure 12: (a) shows the plot of L (normalized) as a function of K with r? = 2 for different 
values of a when the dipole is parallel to its velocity and both lie along the anisotropic direction, 
(b) shows the plot of properly normalized V as a function of L with rj = 2 for the same set of 
values of a. 



(and, in fact, (71) with K replaced by its minimum value K min ). Equipped with this much 
information we can now obtain the plots for the dipole separation and the potential which 
are given in Fig.ll and 12 for 77 = 1 and 77 = 2 (corresponding to v — 0.964) respectively. 
The plots are very similar to those in Sec. 3.4 and so we refrain from giving a detailed 
description. However, note that now the effect of anisotropy is made more conspicuous 
by the significant deviation of the curves from the corresponding isotropic ones. As in the 
previous section, here too, we observe the appearance of a minimal value of the dipole 
separation for the upper unstable branch arising out of the lower bound that was clamped 
upon K. 



4 Comparison among the different cases 

In the previous section we have computed the Q-Q separation and the Q-Q potential 
for different orientations of the dipole and its velocity. Before concluding, let us do a 
comparative study of the effects of anisotropy in all the cases. In Fig. 13 we have given 
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Figure 13: (a) shows the plot of L (normalized) as a function of K with 77 = 0, a = 0.6 for 
different orientations of the dipole and its direction of motion, (b) shows the plot of properly 
normalized V as a function of L with r/ = 0, a = 0.6 for the same set of orientations and direction 
of motion. 



the L(K)-K and the V(L)-L plots for the three surviving cases 8 for 77 = and a = 0.6. 
The legend in the figure needs a little explanation. While the blue line indicates the 
isotropic curve, 'perp' indicates the dipole is moving in the transverse plane, perpendic- 
ular to the direction of anisotropy and 'para' denotes the case where the dipole presents 
itself along the anisotropic direction. While the presence of anisotropy makes itself felt 
in both the cases, the dipole is more affected when it is aligned parallel to the direction 
of anisotropy so that one can write, L max (para) < L max (perp) < L max {isotropic) and 
Vagotropic < Vperp < Vpara- This observation corroborates the findings in [44]. In Fig. 14 we 
have plotted the same quantities, now for 77 = 2 and for all the configurations considered. 
Before delving into the details of the plots, let us again clarify the legend used. Note that 
now there are two isotropic plots, denoted by 'perp' and 'para' indicating the cases where 
the dipole lies perpendicular and parallel to the direction of motion respectively, (ij) 
denotes the configuration where the dipole moves along x % and is aligned along xK Basi- 
cally, one can distinguish between two sectors: one in which the dipole is perpendicular 

8 A little deliberation shows that in the static limit many of the cases collapse into each other and 
need not be considered separately. 



30 



to its velocity (this contains 'perp', (12), (13), (31)) and the one where it is parallel to its 
velocity (comprising of 'para', (11), (33)). The general observation is that the screening 
length diminishes and the potential is weaker for all the cases (ij) shown compared to 
the corresponding isotropic cases. The cases (12) and (31) which merged in the static 
case now splits up and we find that (31) is severely affected when the combined effects 
of velocity and anisotropy are taken into account. This is evident both from the L(K)-K 
and the V(L)-L plots. For this configuration L max drops drastically and also the rise in 
V(L) is appreciable. As discussed earlier too, this is accounted for by the presence of the 
0(y 2 ) term in the anisotropic contribution, which makes the effect of anisotropy quite 
pronounced in this configuration. Both the (13) and (12) cases are mildly affected when 
effects of velocity and anisotropy act in conjunction. For these cases L is slightly sup- 
pressed from the isotropic value whereas V(L) registers a small increase. Turning to the 
other sector, we see that in (11), L decreases marginally which is accompanied by a corre- 
sponding small increase in the interaction potential when we introduce anisotropy and the 
velocity parameter together. However, the (33) plots show a significant departure from 
the isotropic case. For the unstable high energy branch of the potential, the minimum 
allowed separation L min decreases in the order, L min (33) < L min {ll) < L rnin {isotropic). 
On the whole, the plots suggest that the dipole separation and the potential are affected 
the most when the dipole moves along the anisotropic direction (both for the perpendic- 
ular and the parallel orientation and we hope, it will hold true for an other orientation in 
between these two extreme cases), and irrespective of the configuration, the presence of 
anisotropy makes the dipole more susceptible to dissociation. 

In this context, it will be interesting to compare our observations with those extracted 
from other models of anisotropic plasma. In [59] the heavy quark-antiquark static poten- 
tial was computed in an anisotropic plasma employing the hard thermal loop approach. It 
was found out that the presence of anisotropy reduces the screening so that the potential, 
in general gets strengthened and approaches the vacuum potential. The deviation from 
the isotropic potential increases as the value of the anisotropy parameter is increased. Fur- 
ther, the effects are strongest when the dipole is aligned along the direction of anisotropy 
(this is in agreement with our findings here). However, the direction of the shift in po- 
tential (i.e. decrease or increase from the isotropic case) are opposite here (and in [44]) 
and what was found in [59]. By introducing the velocity, we have shown here that for 
sufficiently large value of the velocity the effects of anisotropy on the dipole moving along 
the anisotropic direction will be the strongest. It might be interesting to attempt a similar 
study in the perturbative approach and see if the introduction of the dipole velocity leads 
to results similar to that presented here. Another difference is that the results of [59] 
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Figure 14: (a) shows the plot of L (normalized) as a function of K with rj = 2, a = 0.6 for 
different orientations of the dipole and its direction of motion, (b) shows the plot of properly 
normalized V as a function of L with r/ = 2, a = 0.6 for the same set of orientations and direction 
of motion. 

hold for length scales ~ = l/run where mo is the Debye mass. At this length scale 
the Coulomb part of the potential dominates over the linear confining part. However, 
our analysis here is not constrained in this aspect and in fact, in the static case we have 
identified a range L p < L < L max where the confining part dominates over the Coulomb 
one. Of course, it will be naive on our part to read too much into these comparisons, since 
the physical models in the two cases are completely different, the primary difference being 
that our analysis holds for strongly coupled theories whereas the perturbative calculations 
are valid only in the weak coupling limit. Another curious comparison can be drawn with 
the results of [40] where, among other things, the velocity-dependent Q-Q separation 
and the Q-Q potential were found out in thermal, non-commutative Yang-Mills plasma 
(NCYM). The presence of non-commutativity breaks the isotropy and this is reflected in 
the background metric, where the x 1 direction is taken to be the anisotropic one. Thus it 
is tempting to explore whether the NCYM can serve as a toy model to study the effects 
of anisotropy. The configuration considered in [39] corresponds to that in Sec. 3. 3 in the 
present paper. The counterpart of the anisotropy parameter is the non-commutativity pa- 
rameter in [39]. Surprisingly, we find that the gross characteristics of the results are the 
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same in both models. In NCYM, too, increasing the non-commutativity parameter leads 
to a depreciation of the screening length and a weakening of the interaction potential. Of 
course, there were some additional features in NCYM (like the presence of a cut-off) not 
encountered here. 

5 Conclusion 

Finally, we conclude with a brief summary of the results obtained in the paper. We have 
found out the velocity- dependent Q-Q separation and the Q-Q potential in a strongly cou- 
pled anisotropic plasma at finite temperature via the gauge/gravity duality. The gauge 
theory we take is a deformation of the M = 4 SYM and we take the gravity dual as 
proposed in [41,42]. Barring the screening length in a special case, in all the other cases 
we presented numerical results. The general observation is that when we turn on a small 
value of the anisotropy parameter, the screening length (L max ) decreases and the Q-Q 
interaction becomes weaker so that the dipole becomes more prone to dissociation. We 
considered five different cases, depending upon the direction of motion of the dipole and 
the direction along which it is aligned. While the generic features of the plots are the 
same in all the cases, the minute details vary from case to case. In particular, when 
the dipole lies along the direction of anisotropy the effects are manifested more promi- 
nently in the static case. However, for finite velocity, it is the dipole moving along the 
anisotropic direction that is affected the most. We also set the rapidity parameter rj = 
and recover the static Q-Q separation and the static Q-Q potential. In these cases, our 
findings are consistent with those recently obtained in [44]. Finally, we also compared 
the results obtained in this model vis-a-vis some other models. In particular, we found 
that the results for the static dipole potential provided here (and also in [44]) are different 
from those obtained using standard perturbative field-theoretic techniques in the weakly 
coupled regime. 
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